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SUMMARY 


This report deals with the statistical characteristics 
of surface profiles measured with a stylus tracer; there 
definitions, an application and enhancement using software 
to compensate for stylus geometry effects. After defining 
some of the common height sensitive profile statistics, they 
are used classify the journal surfaces of diesel engine crank 
shafts produce by manufacturing methods that yield signifi- 
cantly different service life. 

Software methods are presented to try to reconstruct a 
surface profile from discrete measurements by accounting 
for the finite radius of the stylus tracer. 

Results indicate that using three parameters: RMS 

roughness, skewness and kurtosis, and a classification 
method termed "separated subspaces", the journal surfaces 
produced by different combinations of grinding and lapping 
can be classified. The work on compensating for stylus 
geometry, which is verified using both mathematical 
simulation and experimental measurements, indicates that, 
at least for simple profile geometries compensation for 
stylus radius' can reduce errors to less than .4%. 
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1. INTRODUCTION 


Before 1950 surfaces were generated and refined exclusively 
by traditional mechanical methods (e.g. cutting, grinding, honing, 
lapping and polishing) . The need to machine the high tempera- 
ture alloys used in ]et engines stimulated the development of 
many additional material removal techniques. These non traditional 
methods include electrochemical cutting and grinding, spark 
discharge machining, electron beam, ion beam, laser and plasma arc 
removal techniques [1]. Thus in the early stages, finish was speci- 
fied by the process, but with these non-traditional methods as the 
drj-ving force, more quantitative methods of specification were needed. 

At the submicroscopic level most surfaces are far from 
smooth and plane, they have the characteristics of a range of 
mountains with peaks and valleys. A number of causes contribute 
to the roughness. First is the mark left by the tool or grit 
itself, which will be of a periodic nature for cutting process 
and more random for abrasive or nontraditional processing 
methods. Second there is a finer structure due to tearing of 
the metal during machining, the debris of the built-up edge and 
the small irregularities in the shape at the tip of the tool. 

Finally, especially in alloy steels there may be microscopic 
cracks at gram boundaries [2] . Thus the resulting surface is 
a function of the process used, the conditions at the cutting 
edge and the material being processed. 

These characteristics of surface roughness are very impor- 
tant in many respects from both a scientific and industrial 
point of view . Particularly in contact problems that 
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involve wear, lubrication, heat transfer and sealing, surface rough- 
ness plays a role [3-12]. For example, in forced convection heat 
transfer, it is well known that the heat transfer rate can be 
enhanced through proper changes of surface roughness [8] . Also the 
rate of fouling or the deposition of scale on the surfaces affects 
the useful life of the heat transfer equipment significantly 
[7]. A large number of engineering components and devices are 
directly dependent upon surface characteristics for their per- 
formance. These include both sliding and rolling bearings of 
all types, seals, brakes, clutches, joints, springs, fasteners, 
cams, splines and gears, particularly if the requirements of 
interchangeability of machine elements considering the fits, 
wear, lubrication, etc., that are involved [13]. 

For some of these applications there is an optimum surface. 

For instance, the cylinder walls of an internal combustion 
engine may be too smooth to allow rapid spreading and wetting 
by oil or too rough to enable the surface asperities to support 
the applied loads without galling [13] . The topic of quantita- 
tively expressing the extent of roughness of the surface is 
really worthy of careful study. Specifically, it would be 
desirable to characterize the form of a surface, be able to 
quantitatively relate this form to the function of the surface, 
and then to know exactly what processes can be used to generate 
this form. 

The ways of measurement of the surface profiles can be 
categorized as follows: 

a. Non-Contact Profile Measurement Methods [14]: 

One example is optical methods, which allow the 
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specimen to be investigated without destroying it or 
subjecting it to a strain or wear. Other methods that 
have been used involve capacitance or pneumatics as the 
measurement principle, 
b. Contact Profile Measurement Methods; 

Instruments with stylus tracers in mechanical contact 
with the surface are the most common way of measuring surface topo- 
graphy. Stylus methods have a shortcoming in that there 
exists an error from the influence of stylus geometry, 
but it does provide an immediate numerical characteriza- 
tion of a surface, so it is used widely in industry, is 
the most direct measurement of geometry and is used in 
U.S. Standards [15]. 

Statistical considerations are intimately tied up with the 
measurment of surfaces. Statistical parameters are used to 
characterize different surfaces with the expectation that there 
will be little variation in these parameters over the surface 
[16]. Various modifications and improved surface finish para- 
meters have been proposed by Reason [17] , pesante [18] , Ehrenreich 
[19], Teague [20] and others. 

This report concentrates on measurements made with stylus 
devices that are digitized representations of a stylus trace. 

As a starting point, some of the common parameters used to 
characterize the form of a profile trace are defined. This is 
followed by an application of these parameters to the problem 
of characterizing the journal surfaces of crank shafts that are 
produced by different manufacturing methods and have vastly 
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different life in service. The final piece of work develops 
new methods for deconvolving or compensating for the effects of 
stylus geometry on the measurements made with a tracer. 
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II. STATISTICAL CONCEPTS FOR DESCRIBING TYPICAL SURFACE 


TOPOGRAPHY PARAMETERS 

The most common surface finish measurement variable is 
roughness height, which is the numerical value of the average 
distance, in micro-inches or microns, of each point on the 
surface profile from a defined line called the reference center 
line. Once this reference center line is set, each roughness 
height measurement, Y^ for i = 1,...,N, of the surface are 
referenced normal to this line. All computations of the charac- 
teristics of the profile are based on the measured roughness 
height. Obviously it is essential to define a reference 
center line properly. 

There are several methods that have been used in defining 
the reference center line [21,22], they are: 

a. Envelope Method (E - System) : 

Imagine that there is a large circle (often 25 mm 
in diameter) rolling over a surface, and regard 
the locus of the center of this circle as the reference 
center line. This method is used in some European 
countries. 

b. Ten Point Average Method: 

This method requires finding the 5 highest peaks and 
5 lowest valleys of the profile, and calculating the average 
value of these ten points. This average value is regarded 
as the reference center line. 
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c. Mean Line Method; 

The mean line is selected so that on each side of it the 
areas enclosed by the profile are equal, i.e. the centroid 
of the profile. For discrete profiles, the area for 
each profile is assumed to be a rectangle of height 
and with a constant width Ax. It turns out that we can 
use an alternative definition for easy computation which 
takes the reference center line as a line parallel to 
the general direction of the profile such that the average 
height of the profile on one side of it is equal to the 
average height on the other side. With the previous 
assumptions, the mean line in this case is simply the 
average, i.e. if the height of the point at is Y^, 
the mean can be mathematically expressed as : 


Y = 


1 

N 


N 

E 

i=l 


Y. 

1 


( 1 ) 


This method is the standard in U.S., Canada and Britain, 
d. Least Squares Line Method: 

The well known formulas for linear regression are 
used to get the least square line, which is regarded as 
the reference center line. With this method the refer- 
ence line IS a function of position as given by: 

Y^ = a + bx^ (2) 


where 

a = the intercept of the least square line with 
the y axis, and 

b = the slope of the least square line (3) 
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In practice, the Envelope Method is not used very often 
because of the difficulty of detemining the locus. The Ten 
Point Average Line Method leads to a reference center line 
below the major surface features for deeply pitted surfaces^ 

As a result it is common that people often use the Mean Line 
Method and the Least Square Line Method. 

The Mean Line Method gives a "horizontal" reference line 
which cannot compensate for the "tilt" in the profile, whereas 
the Least Square Line Method does compensate for tilt in the 
experimental setup. 

Now, let us have a brief survey of the statistics which 
are most commonly used to represent the properties of a 
measured surface. All of these parameters are based on a 
profile like that shown in Fig. (la) . They refer to deviations 
from a reference line based on one of the methods previously 
described. Therefore, all computations are made using; 

y = Y - Y (4) 

where ? is the reference line, Y^ is a measured value and 
y IS the deviation from the reference. This transformation 
leads to another discrete profile which may be interpreted as 
shifting the measured profile to a zero mean level, as shown in 
Fig. ( lb) . 

HEIGHT PARAMETERS: 

Measures of dispersion show the degree of spread of 
the data around the central value. The most common one 
IS standard deviation, or RMS roughness. Based on devia- 
tions from the mean given by Eq. (4) , the RMS roughness 
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IS defined as: 


R 

q 



2 


(5) 


Today, because of its greater simplicity, arithmetic 
averaging is much more commonly used, and is, in fact, 
the American Standard for roughness. Arithmetic average 
(AA) roughness is defined as: 


R 


a 


1 '' 
— z 

N ^ 


i=l 


HEIGHT DISTRIBUTION PARAMETERS: 


(6) 


There are two parameters which were proposed by Al- 
Salihi [23] to describe height distribution, they are 
skewness and kurtosis. While they are used in charac- 
terization of surface profiles, they are well known 
descriptors of statistical distributions. 

Skewness means lack of symmetry, and measures of 
skewness show the extent to which the distribution departs 
from symmetry. Skewness is defined as: 


Yl 


1 

N 


N 

Z 

1=1 



3 


(7) 


Refer to Fig. (Ic) . If = 0, the distribution is 
symmetric, such as a Gaussian distribution, shown as 
curve 1. If > 0/ the distribution is skewed to a 
higher level as shown by curve 2. Whereas, if < 0, 
the distribution is skewed to a lower level as shown in 
curve 3. The positive skewed surfaces (y^ > 0) is thought 
to be more suitable for load carrying than surfaces 
negatively skewed. 
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Kurtosis may be defined as "peakness". A measure 
of kurtosis serves to differentiate between a flat 
distribution curve and a sharply peaked curve. In 
other words, it enables the squareness of the profile 
to be described. Kurtosis is defined as: 


^2 


1 

N 


N 

S ( 
1=1 



4 


( 8 ) 


For a Gaussian distribution, Y 2 equal to 3, which 
IS shown as curve 1 in Fig. (Id). If Y 2 ^ distri- 

bution is more sharply peaked than Gaussian as shown in 
curve 2, and is defined as leptokurtic. If Y 2 
distribution is flatter than Gaussian as shown in curve 
3 and is defined as platykurtic. 

LENGTH SENSITIVE PARAMETERS : 


One parameter in this group is the autocorrelation 

function, which was first noticed by Wormersley and 

Hopkins [24] as a time series. However it was Peklenik 

[25] who first applied it to classification. 

The autocorrelation gives an estimate of the relation 

between y and y , , which are the values of y. at hori- 
-'i -'i-k' ■'1 

zontal intervals of length, k(Ax). Autocorrelation is 
defined as: 


A. 



N 

E 

i=k+l 

N 

E 

i=l 


^i^i-k 


Y. 


( 9 ) 


in addition to the autocorrelation, the spectrum. 
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whxch is the Fourier transform of the autocorrelation, is 
often used. Often the spectrum is most effective when 
dealing with highly periodic profiles. 
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III. HOW ItSASURSMENTS APE MADS 

The procedure used to make profile measurements involves 
digitizing the analog stylus deflections, and storing this 
data. The setup used is typical of many laboratory installations, 
this one using a Bendix Proficorder and a Digital Equipment Corpor- 
ation (DEC) LSI-11/2 microcomputer. 

Figure (2) shows the stylus transducer setup that is used 
to convert the vertical motion of the diamond stylus, with 
radius r, to an analog voltage. The lever arrangement causes 
the core of a linear variable differential transformer (LVDT) , 
to move. The resulting A.C. signal is demodulated to provide a 
D.C. voltage proportional to the deflection of the stylus tip. 

The stylus traverses at a linear velocity V, which is 
assumed to be constant. To assure a straight path for this 
motion, the stylus is referenced to an optical flat. However, 
this straight path does not assure that the stylus has a path 
parallel to the surface being measured, so that it is possible 
to have a "tilt" in the measured profile. 

The data acquisition setup is shown in Fig. (3) , and is 
designed to provide an analog trace of the surface, as well as 
a digitized trace. A Brush recorder is used to indicate the 
analog trace on one channel, with the signal coming directly 
from the tracer amplifier. Between the amplifier and the analog 
to digital converter (ADC) on the microcomputer, an active low 
pass filter is installed to avoid aliasing as explained below. 

The aliasing problem can best be explained in terms of the 
sampling interval. Ax. If we sample at points which are too 
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Figure 2. Schematically illustrating the (a) linear variable 
differential transformer, (b) locus of stylus tip when measuring 






A 

B 

C 


Journal Specimen 
Stylus 

Linear Variable 

Differential 

Transformer 


D — Amplifier 
E — Filter 
F — Recorder 

G — LSI- 11/2 Microcomputer 
H — CRT 


Figure 3. Schematically illustrating the 
set-up. 


15 







close together, it will yield correlated and highly redundant 
data, and thus unnecessarily increase the labor and cost of 
calculations. On the other hand, sampling at points which are 
too far apart will lead to confusion between the low and high 
frequency components in the original data. This later problem 
is called aliasing. 

Consider a continuous record which is uniformly sampled 
with h seconds time interval, i.e. a sampling rate of 1/h 
samples per second. If the velocity of stylus motion is V, the 
sampling interval will be Ax = V*h, refer to Fig. (2). However, 
we need at least two samples per cycle to define a frequency 
component in the original data. Hence, the highest frequency 
which can be defined by sampling at a rate of 1/h samples per 
second is l/2h Hz. Frequencies in the original data above 
l/2h Hz will be folded back into the frequency range from 0 to 
l/2h Hz, and be confused with data in this lower range. This 
cutoff frequency is called the Nyquist frequency. To be on the 
safe side the filter break frequency is set at l/3h. 

The digitized signal is sent back to the second channel of 
the Brush recorder using a digital to analog converter (DAC) so 
that it can be visually compared with the incoming analog 
signal. For subsequent analysis, the digitized data is also 
stored on a floppy disk. 

Once the data on a profile has been acquired, it is condi- 
tioned as follows. The trend is defined as any frequency com- 
ponent whose period is longer than the record length. This 
type of component cannot be removed by highpass digital filtering 
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as will be mentioned later. Here we chose the least squares 
procedures, Eqs. (2) and (3), to remove the linear trend, which 
usually arises from "tilt" or lack of parallelism between the 
optical flat and surface being measured. 

To remove waviness often associated with errors of form, 
highpass filtering, like that often done with wavelength cutoff 
analog circuitry, is used. This can be done by fast Fourier 
transform because only a finite range Fourier series or transform 
can actually be computed with digitized data, and this finite 
range can always be considered as the period of an associated 
Fourier series. Digital filtering methods are used to filter 
out the lower frequencies, (long wavelength waviness) of the 
profile by choosing a wavelength cutoff. 

Once this is done, the profile statistics described in 
Section II can be computed. 
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IV. AN APPLICATION OF SURFACE STATISTICS TO CLASSIFICATION 


OF JOURNAL SURFACES ON CRANK SHAFTS 

The normal finishing steps on journal surfaces involves 
grinding and lapping. The relative direction of these two 
operations is felt to be critical. For example, if the crank 
will rotate clockwise, then the grinding should be done counter- 
clockwise followed by lapping in a clockwise direction. One can 
speculate that this order could tend to minimize the directional 
tendency of asperity tips, i.e. the grinding may give the asperi- 
ties a direction and if lapping works on the tip of the asperity 
It will flatten the asperity and shift material in the opposite 
direction. It is said that there is quite a difference in bear- 
ing life v/hen using the journal bearings made through different 
manufacturing methods. For example, the life of the journals 
which are ground and lapped in the same direction is 500 working 
hours, whereas the life of those which are ground and lapped in 
the opposite direction is about 5000 working hours. The effect 
of the finishing steps is so great as to be worthy of studying. 

An experiment was made by applying the previous ideas on 
the measurement of surface profiles to several journal surfaces 
of diesel engine crank shafts. An additional purpose was to 
determine if there is a parameter or several parameters that 
can be used to classify the journal surfaces according to their 
manufacturing steps and the relative directions of grinding and 
lapping. 
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SAMPLE SPECIMENS: 


Coupons to be measured were cut out of the journal 
surfaces of new engine crankshafts, some from the main 
bearing surfaces and others from the connecting rod 
throw. Figure (4) , illustrates where the various samples 
are located on a typical crank originally. 

Three sets of specimens, from three cranks, were 
made available; ground only and unlapped (U), ground and 
lapped in the same direction as grinding (LSD) , and ground 
and lapoed in the direction opposite to grinding (LOD) . Table 1 
groups the samples according to their manufacturing procedure. 
MEASURING : 

The set up of this experiment is the same as shown 
schematically in Fig. (3). Surface profile traces 
were made using a Bendix Proficorder equipped with a 
stylus having a 12.7 ym radius. The analog output of 
the stylus displacement was digitized, bypassing the 
analog filters used for setting the wavelength cutoff. 

A Krohn-Hite 3323 active filter acted as an anti-aliasing 
filter. Based on the traverse speed of the stylus, .3175 mn/s, 
and selection of a spatial sample interval, AX = .005 mm, 
the temporal rate is determined and the break frequency 
for the antialiasing filter was selected on the conservative 
side to be one third the sampling frequency. 

Using the conditions mentioned above, 4 longitudinal 
traces were made at different positions on each journal 
coupon. In each case, 512 points were sampled in each 
trace, for a total stroke of slightly more than 2.5 mm. 
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Rod #3 
( E4,E2 ) 


Mam #5 
( C4,D4,G4 ) 
( C3,C],G] ) 



Rod #4 
( B4 ,f4 ) 
( F3,F2 ) 


Figure 4. Schematically illustrating the 
location of the samples on a crankshaft. 
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hanuf acturind 
Procedure 

Specimen 

Identification 

Ground only 
Unlapped < U ) 

B4»D4fE:4»E2 

Ground and 
Lapped in Same 
Direction (LSD) 

C4»C3»C1 

Ground and 
Lapped in Opposite 
Direction (LOD) 

F4»G4»F3rGl»F2 

r 


Table 1. Group the samples sccordinsi to their 
manuf acturins procedures 
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These data were analyzed usxng the least squares 
reference line, Eq. (2) , and the statistical parameters 
given by Eqs. (5) through (8) were computed. Reference [27] 
g,iv,es complete details on the computational methods used. 

THE IDEA OF THE "RMS - SKE^VNESS - KURTOSIS SPACE"; 

By comparing all the statistical parameters listed 
in Tables (2.1-2. 3), it is true that the "ground only" 
specimen can be easily recognized from the other two 
kinds of specimen by just looking at the arithmetic 
averaging or RMS roughness. Owing to the similarity 
between arithmetic average and RMS, we 
choose only one of them, Rq, as a charac- 
teristic. The remaining parameters to describe the 
characteristics of the profile are Rq, 
skewness and kurtosis. Because every specimen has a 
set of values, we can regard it is a set of coordinates 
in a space constituted by these three characteristic 
axes. Since the values of coordinates are related to 
the wavelength we choose, we hope we can classify all the 
data points located in the defined space into three groups 
by choosing a suitable wavelength cutoff. 

THREE MODELS OF CLASSIFICATION: 

It is supposed that all the data points with the same 
manufacturing procedure will cluster into a sphere around 
a certain center. We took the average of all the data 
points with the same manufacturing method as the center of 
the sphere, shown in Table (3), and found that wavelength 
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Wavelen^ bh 
Cutoff 


0*8 miTi 


Table 2*1 


Group 

Type 

ID 

Ra 

F<io 


•^4 


B4 

0*912 

1*099 

-0*266 

2*836 

U 

D4 

0*579 

0*719 

-0*163 

2*785 

E4 

0*99A 

1*251 

0*014 

3.008 


E2 

1*387 

1 ,745 

0*412 

3*249 


C4 

0*218 

0.266 

0*126 

2*934 

LSD 

C3 

A , nnn 

0.301 

-0*828 

8* 132 


C1 

0*283 

0*365 

-0*053 

5*273 


F4 

0*387 

0*457 

0*009 

2*160 


G4 

0*229 

0*275 

0*081 

2*431 

LOD 

F3 

0*404 

0*499 

-0*334 

3*420 


G1 

0*330 

0*447 

-0*610 

4*359 


F2 

0,411 

0*538 

0,-445 

4*008 




Statistical parameters for 0»8 um wavelength cutoff 
(each parameters Sot from the averase of 4 traces) 



f 


Wavelength 

Cutoff 


0«25 miTi 


Table 2.2 


Group 

Type 

ID 

Ra 

Ro 




B4 

0.579 

0.755 

-0.485 

4.204 

U 

D4 

0.561 

0.706 

-0.234 

3.044 

E4 

0.877 

1.103 

-0.161 

3.163 


E2 

0.969 

1.218 

0.002 

3.113 


C4 

0.165 

0.211 

-0.140 

3.745 

LSLi 

C3 

0.134 

0.188 

-0.339 

9.768 


Cl 

0.163 

0.216 

0.615 

16.867 


F4 

0.269 

0.338 

-0.033 

3.140 


G4 

0.154 

0.199 

0.260 

3.734 

LOD 

F3 

0.315 

0.405 

-0.245 

3.992 


G1 

0.178 

0.246 

-0.408 

6.222 


F2 

r 

0.278 

0.388 

0.138 

f 

4.893 


Statistical parameters for 0.25 mm wavelength cutoff 
(each parameters sfot from the average of 4 traces) 




Uavel ensth 
Cutoff 

Group 
r ype 

ID 

Ra 

Rq 


I 



B4 

0.441 

0.570 

-0.385 

♦ 

03 

(> 


U 

D4 

E4 

0.451 

0.738 

0.567 

0,934 

-0.169 

-0.002 

2.945 

3.140 



E2 

0.792 

0.996 

-0.096 

3.074 



C4 

0.084 

0.106 

-0.217 

3.233 

0«08 mm 

LSD 

C3 

0.084 

0.118 

-0.377 

8.774 



Cl 

0.096 

0.134 

0.690 

23.461 



F4 

0.135 

0.175 

-0.178 

3.869 



G4 

0.102 

0.133 

-0.125 

4.249 


LOD 

F3 

0.175 

0.233 

-0.854 

5.725 



G1 

0.106 

0.156 

-0.231 

11.198 


f , 

F2 

f i 

0.147 

f _ . i 

0.200 

f . i 

-0.366 

f i 

4.739 

f 


Table 2«3 Statistical parameters for 0*08 mm wavelensth cutoff 
(each parameters sot from the average of 4 traces) 



f 


Wavelength 
Cutoff 
( mm ) 

U 

LSD 

LOD 

0,8 

(1.204»-0.001» 2.970) 

(0.311»-0.252» 5.446) 

(0.443»-0.082» 3.277) 

0.25 

(0.946r-0.220» 3.381) 

(0.205» 0.045»10.127) 

(0.315r 0.058» 4.396) 

0.08 

(0.767»-0.163» 3.3A1) 

(0.119» 0.032»11.823) 

(0.179f-0.351r 8.028) 




Table 3 Reference center points (Standard deviationr Gkewness» Kurtosis) 
(each value is the average of the same aroup) 



cutoff did affect the position of this center. Thus we 
examined some models to find the best wavelength cutoff 

for clearly distinguishing the three conditions, 
a. Totally Separate Sphere Range Model: 

If we choose the distance between the farthest 
individual point and the corresponding center point 
as the radius and draw a sphere , we get three 
spheres with three different centers. The optimal 
condition for which we can distinguish these three 
spheres, which stand for three different ways to 
make journal bearing, is to maximize the distance 
between all the centers. 

Since we have three centers, the distance between 
every pair of centers are listed in Table (4) , the 
maximum radius each sphere may have can be considered 
as follows. Say we have three spheres with centers 
at points A,B,C, and the corresponding sides are a,b,c. 

If there is a smallest side, e.g. c, then both spheres, 
which have their center at either tip of side c, may have 
a maximum radius equal to c/2. The maximum radius of 
the third sphere is equal to the difference between the 
smaller side and c/2. This can be easily understood 
when v/e look at the triangle constituted of the three 
centers as shown in Fig. (5). 

Following the previous idea, the procedure to do 
this IS to find the three centers under different 
wavelength cutoffs. Then calculate the maximum radius 
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Wavelensfth cutoff (mm) I 

. 1 



0.8 ! 0.25 

.1 

1 

1 

1 

! 

CO 1 
O 1 

♦ 1 

O 1 
1 
1 
1 
1 
1 

u - 

LSD 

9.599 : 10.361 

. 1 

10.831 

LSD - 

LOD 

3.058 ! 5.837 

.... 

5.424 

LOD - 

U 

• • 

7.652 J 6.960 

7.733 


/_ / ^ 


Table 4 Distance between sroup centers 
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each sphere may have, and check the number of data 
points that fail to fall into the corresponding 
range. The wavelength cutoff which minimizes this 
failure is the one selected. 

b. Minimized Overlapping Area Model: 

This idea is somewhat similar to the previous 
one. The main difference is that we choose the 
longest distance between each data point and its 
corresponding center as radius of the sphere. We 
got three spheres from three cases. The wavelength 
cutoff we need is the one which produces minimum 
overlap of the three spheres, as shown in Fig. (6). 

c. Separated Subspaces Model: 

The third approach is based on the idea that if 
we find the three centers first, then the data points 
of the same group should have a shorter distance from 
the corresponding group center than those from the 
other two group centers. This can be expressed geome- 
trically, refer to Fig. (7). Imagine a triangle with 
three centers A,B,C as its tips. The three planes 
which are perpendicular and bisect the three 
sides individually will intercept at a line called 
the centroid line. These three planes divide the 
space into three subspaces. The data points from the 
same group should fall into the same subspace. Since 
the distance between a data point and the center point 
in the same subspace will be the shortest one among 
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the three possible alternatives. The classification criteria is 
to select the region which minimizes the distance to the corres- 
ponding center. The best wavelength cutoff is the one for which 
the most data fit the model. Table (5) shows the results based 
on the totally separated sphere criteria, while Fig. (6) graphic- 
ally illustrates the minimized overlapping area idea. The re- 
sults in Table (5) suggest that the .25 mm cutoff gives the best 
classification because the number of correct classifications is 
greatest. Figure (6) also suggests that the .25 mm cutoff is 
the best to use, because the overlap area is the smallest. 

Results with the third model, that using the separated sub- 
spaces idea, are given in Tables (6. 1-6. 3). The distances to all 
three center points for each specimen are given in each row, 
with the selection based on the shortest distance. The last 
columns in Tables (6. 1-6. 3) indicate a correct or incorrect 
classification. Based on these results, again the .25 mm cutoff 
has the greatest discriminating power. 

We may conclude that among the three models mentioned 
above, separated subspaces model is the most suitable one to 
classify these journal surfaces with regard to their manufactur- 
ing method. Also, 0.25 mm wavelength cutoff is proved to have 
a better power to subdivide surface roughness of journal surfaces. 
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1 

» 

Within or out of the predicted ransie ! 

at different wavelensth cutoff <mm) ! 

Type 

ID 

4 

0.8 

0.25 

> » 
0.08 


B4 

IN 

IN 

IN 

U 

D4 

EA 

IN 

IN 

IN 

IN 

IN 

IN 


E2 

OUT 

IN 

IN 


CA 

OUT 

OUT 

OUT 

LSD 

C3 

our 

OUT 

OUT 


Cl 

OUT 

OUT 

OUT 


F4 

IN 

IN 

OUT 


G4 

OUT 

IN 

OUT 

LOD 

F3 

OUT 

OUT 

OUT 


G1 

OUT 

OUT 

OUT 


F2 

OUT 

IN 

OUT 


— — — - ' / - 


Table 5 Results of apelyinS the idea of totally separated 
sphere ranse 
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Type 
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LOD 


Preh icteri 
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Correct ( C ) 
or 

Wron'3 (W) 



D4 

2 . 854 

8.306 

6.026 

1) 

C 


D4 

5.113 

4.954 

2.917 

LOD 

w 

u 

E4 

0.499 

10.071 

8.139 

U 

c 


E2 

6.814 

15.956 

13.923 

U 

c 


C4 

9.460 

4.558 

2.752 

LOD 

u 

LSD 

C3 

13.287 

6.359 

9.015 

LSD 

c 


Cl 

8.711 

2.067 

2.163 

LSD 

c 


F4 

7.509 

4.443 

1 .446 

LOD 

c 


04 

9.337 

4 . 504 

2.489 

LOD 

c 

LOD 

F3 

7.807 

2.881 

2.587 

LOD 

c 


01 

9.812 

3.985 

5.392 

LSD 

u 


F2 

8.077 

7.469 

5.402 

LOD 
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Table 6.1 Results of ■Separated Subspaces' with 0.8 mm wavelength cutoff 
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Wave 
Length 
Cutoff 
< mm ) 

Group 

Type 

U 

ID 

D4 

D4 

E4 

E2 

Distanc 
ref erer 

U 

3.370 
2.423 
1 . 694 
3.522 

:es to the i 
ice center f 

LSD 

9.667 

9.114 

1 1 .550 

12.329 

.hree 

'Oints 

LOD 

6.987 

5.060 

8 . 268 

9.136 

Predicted 

Type 

ti 

LI 

U 

U 

Correct ( C ) 
or 

Wron^ (W) 
C 
C 
C 
C 



C4 

7.397 

6.646 

2.327 

LOD 

W 

0.25 

LSD 

C3 

9.980 

3.864 

6.798 

LSD 

C 



Cl 

17.457 

8.826 

13.696 

LSD 

C 



P4 

6.359 

7.155 

1.565 

LOD 

C 



G4 

8.879 

6.774 

2.426 

LOD 

c 


LDD 

F3 

5.445 

7.076 

3.182 

1 OD 

c 



G1 

7 . 782 

5.997 

5.049 

LOD 

c 



F2 

6.793 

5.622 

1.193 

LOD 

c 




Table 6.2 Results of "Separated Subspaces" with 0.25 mm wavelength cutoff 
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U 

U 

U 
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or 
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C 
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LSD 

C3 

8.716 

5.101 

1.001 

LOD 
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Cl 

22.733 

13.370 

18.620 

LSD 
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F4 

5.941 

8.245 

4.504 

LOD 

c 


G4 

6.411 

7.736 

4.427 

LOD 

c 

LOD 

F3 

9.046 

10.815 

5 . 560 

LOD 
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01 

9 . 959 
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3.397 
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Ul 


F2 

6.176 

8.165 

3.299 

LOD 
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Table 6«3 Results of "Separated Subspaces* with 0.08 mm wavelensth cutoff 



V. THE COMPENSATION OF MEASURED SURFACE PROFILE 


If we look at the Fig. (2) , we may see that when the 

stylus moves on the profile, the height we really measured at 

position IS the locus of the stylus center, 

which IS Y Y.' distant from the real height Y.. Because of 
11 1 

this inevitable measuring error inherent from the geometry of 
the stylus, particularly the finite radius r, the profile we 
measured is only the locus of the stylus center, which is 
different from the true profile, shown in Fig. (2). It is for 
this reason that it is necessary to compensate for this error so 
that an actual profile may be drawn. 

In general, the best that can be done is to approximately 
reconstruct the true profile. The following models are those 
we chose to compensate for some of the error. 

As a standard for comparison, the proposed compensation 
or deconvolution methods are evaluated in terms of their 
effects on the height sensitive parameters given by Eqs. (5-8) 
for both mathematically simulated surfaces with known para- 
meters and a measured triangle shaped calibration surface. 
STRAIGHT LINE PROFILE MODEL: 

Imagine that we have an oblique profile PP inclined 
at angle 0 as shown in Fig. (O.a) , and consider 
the tip of the stylus as a ball with radius r. When the 
stylus measures the oblique line in a direction which is 

parallel to the datum line, at position the contact 
point of stylus and surface is C^, the center of stylus tip 
is O . The measured height at position X is Y 
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whoroas thG true height of the profile at oosition 

IS Y , as shown in Fig. (8a), so there exists an error Y.Y ' 

1 11 

between the measured height and the true height. If we 

know the slope of the oblique line and the radius of the 

stylus tip, we can get the actual height at some position 

by subtracting a distance Y^Y^' from the measured height 

at the same position. The distance 

Y Y . ' = r ( sec 0-1) 

1 1 

= r (/l + (dy^/dx^)^ - 1 (10) 

where r is the stylus radius, and 0 is the tangential 
angle at contact point C^. 

CONVEX AND CONCAVE PROFILE MODELS: 

Referring to Fig. (8b and 8c) , imagine that we have a convex 
or concave profile with a stylus running over it. The 
contact point of the stylus and profile is C^ , the center of 
stylus tip at position is is the measured height 

and Y^ is the true height at the same position X^. If 
the radius of stylus tip is r, the radius of curvature of the 
profile at position is R^ with center at 0^' and the 
tangential angle 0 at contact point C^ are known, we may 
get the true height Y^ by subtracting the distance Y^Y^' 
from the measured height Y^', which in this instance gives 

/2 2 2 ~ 

Y.Y ' = + /r. - D.sin 0 + D.cos 0 - r (11) 

11—11 — 1 

where the first term is negative and the second positive if 
the profile is convex and the opposite signs apply for a 
concave profile, and 
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n, = r + R . 

1 1 

The restriction on this solution is 


2 

R > 
1 — 


(dy^/dx^) 

The first model, Eq. (10) contains only two variables but 
the other model, Eq. (11), has more variables that must be 
determined, in addition to the restriction given by Eq. (13) , 
which can be violated when the angle becomes large or when R 
becomes smaller than r. Since the first model is simpler to 
follow and above all, with no limitations on application, 
we shall continue our discussion on modifying processes based 
only on Eq. (10) . 

When applying the straight line model to the measured 
profile, we need to know the two variables first, the radius 
of the stylus and the slope of the tangential line at the 
contact point. The former can be measured directly, but the 
slope can only be estimated from the profile measurements. 

To estimate the slope, dy^/dx^, two approaches are used. 
The first, designated Method I, uses a backward difference 


approximation to the derivative 


^1 - yj-1 

Ax 


where the y^ ' s are the measured profile heights cind Ax is the 
sample interval. Method II amounts to a central difference to 
estimate the slope, i.e. 
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2Ax 


( 15 ) 


dy. 

1 


(2) 


^i+l 


1-1 


for approximating the slope. 
VERIFICATION OF THE PROPOSED METHODS: 


To determine how much improvement is possible using 
the approaches outlined in the previous section, both 
mathematical simulation and actual profile measurements 
were used. Fairly simple analytically described functions 
were used for the profile shapes, viz. a sine wave and a 
triangle wave. For these shapes, the parameters given 
by Eqs. (5-8) can be calculated analytically and are: 



Sine 

Triangle 

^a 

H/it 

H/4 

R 

q 

H»/T/4 

h/T/6 

^1 

0 

0 

^2 

1. 5 

1.8 


Where H is the peak to valley height. Note that all these 
pareimeters are independent of the period, meaning they are 
only height sensitive. Furthermore, since the skewness 
and kurtosis are normalized by R , they are dimensionless 
numbers. 

The purpose of the mathematical simulation was to be 
able to eliminate errors introduced in the profile measure- 
ments that can be attributed to the manufacturing of refer- 
ence standards. The analytically defined profile was 
generated, and the simulation program was designed to provide 
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the resulting motion of the stylus as it traversed this 
profile. To do this, suppose that the stylus is now at 
horizontal position X^, and the profile neighboring to this 
position IS decomposed into discrete points. Assume the verti- 
cal position of the stylus center is and calculate all the 

distances between this assumed center point and all the 
discrete profile points. These distances are comoared with 
the stylus radius; these distances 

must be all no less than the stylus radius, with at 
least one distance equal to the stylus radius. We can 
find the measured height by iterating the position of 
the stylus center vertically. 

The simulated triangle shaped profile was based on 
the geometry of a roughness specimen that is used for 
calibration purposes. This standard is certified to 
have a roughness R = 3.124 + .10 ym. A profile was made 
on this standard using the setup shown in Fig. (3) using 
a Ax = .005 mm and a stylus radius of 12.7 ym. Figure 
( 9 ) shows the analog and digitized trace of this sample, 
and this same data was used for subsequent analysis. 

Using the data in Fig. ( 9 ) , the period of the triangle 
wave was found to be P = 93.133 ym, and the peak to valley 
height H = 12.497 ym. 

The simulation used the 12.7 ym radius and the afore- 
mentioned triangle wave characteristics. With the same 
stylus radius, a sine wave with an amplitude and period 
the same as the triangle wave was also simulated. 


43 







Table (7) indicates how the theoretical parameters compare 
with those based on the measurements before compensation to 
account for stylus radius. Table (8) shows how the height sensi- 
tive characteristics are affected by applying Eq. (10) to the 
simulated measurement and using Method I, Eg. (14), and Method 
II, Eq. (15), to estimate the profile slope. The errors that 
correspond to the results in Table (8) are given in Table (9) . 
Table (10) , while similar to Table (8) , differs from it in that 
the sine wave profile is simulated, rather than a triangular 
profile. 

There are obvious differences between those parameters 
calculated from different geometric profiles. The theoretical 
parameters, which are calculated based on the ideal triangular 
profile with the specified height and period, are obviously 
closer to those experimentally measured with stylus. However, 
there is some difference between the theoretical values and 
measurements due to the inability to make a perfect standard. 

The errors in Table (9) indicate that some improvement in 
obtaining the height sensitive parameters can be obtained using 
the methods based on Eqs. (10), (14) and (15). Specifically, 
for triangle shaped surfaces, the corrections based on Eq. (10) 
and Eq. (14) gave the best results, with all errors less than 
.4%. While the overall improvement is not great for an idealized 
surface like a triangle, the improvement could be greater for 
real surfaces that are more random. 
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Ra 

< >im ) 

Rq 

( jjm ) 

K 


3.014 

3.465 

-0.026 

1.767 


/ / y / / 


Table 7« Parameters of measured triangular wave 



Ra 

Rq 

i^. 


Theoretical 

3.124 

3.608 

0 

1 4 B 

Simulated 

Measurement 

3.113 

3.592 

-0.016 

1.776 

Corrected 

usinS 

Eos.(lD) and (14) 

3.118 

3.604 

0.033 

1.793 

Corrected 

ijsins 

Eos. (10) and (15) 

3.111 

3.594 

0.066 

1.788 




Table 8 Comparison for Simulated TrianSle Wave 


46 













y y 



Error 

E rror 

Error 

Error 


of 

of 

of 

of 


Ra 

Ra 

y^ 


Simulated 

~0,011 

-0,016 

-0.016 

-0.024 

Measurement 






" 0*37% 

-0.42% 


1.33% 

Corrected 

0*006 

-0.004 

0 . 033 

-0.007 

usinS 





EIos.(lO) and (14) 

-0 . 20% 

“0.10% 


-0.38% 

Corrected 

-0.013 

“0.004 

0.066 

-0.012 

usind 





Eqs.(IO) and (15) 

-0.49% 

-0.37% 


-0.69% 




Table 9 % Error of parameters with theoretical set as basis 

for simulated brianSular wave 
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f 



Ra 

Ra 



Theoretical 

3.978 

4.418 

0 

1 .5 

Simulated 

Measurement 

3.911 

4.369 

-0.185 

1.547 

Corrected 

usina 

Eqs.(IO) and (14) 

4.038 

4.481 

-0.025 

1.482 

Corrected 

usins 

Fas. < 1 0) and < 15) 

4.028 

4.467 

-0.021 

1.478 




Table 10 


Comparison for simulated sine wave 




CONCLUSIONS 


Work presented in this report was aimed at trying to find 
out which statistical parameters estimated from digitized pro- 
files can discriminate between the surface topography of journal 
surfaces produced by different msuiufacturing sequences. Results 
were also presented on different ways to compensate for the 
errors introduced by having a finite stylus radius on surface 
profile measurements. 

The results presented indicate that: 

1. One parameter or pair of parameters can not sufficiently 
discriminate between different surface topography. Instead, 
using three parameters, i.e. RMS roughness, skewness, kurtosis, 
may do this application well. 

2. The "Separated Subspaces" in RMS - skewness - kurtosis 
space IS a good model to deal with the classification of the 
journal surfaces and may be applied to other categorization 
work. With this criterion, the proper wavelength cutoff for 
classification analysis is 0.25 mm cutoff. 

3. The proposed ways for compensating for the stylus error may 
correct the measured profile and make it closer to the actual 
profile. Using both simulated and measured profiles, it was 
found that improvements can be made using the proposed method, 
particularly when backward differences are used to estimate the 
profile slope. Errors for the height sensitive profile para- 
meters of simulated surfaces were less that .4%. 
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